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Abstract 

We present the application of a new jet reconstruction algorithm that uses a Gaussian filter to 
locate and reconstruct the jet energy to p+p and heavy ion data from the PHENIX detector. This 
algorithm is combined with a fake jet rejection scheme that provides efficient jet reconstruction 
with an acceptable fake rate. We show our first results on the measured jet spectra, and on jet-jet 
angular correlation in p+p and Cu+Cu collisions. 

The study of jet physics at RHIC plays an important role in understanding the parton-medium 
interaction, in-medium fragmentation properties, and applied to the RHIC spin program, also the 
study of the proton structure. We therefore set out to systematically investigate the feasibility of, 
and the appropriate approach of performing jet reconstruction using a limited aperture detector 
like the PHENIX central arms [1] and including heavy ion collision systems. Neither is a well- 
studied aspect of jet reconstruction, and both have been considered as challenging 

We observed that the flat weighting in traditional jet reconstruction algorithms is particularly 
prone to fluctuations at large angle. A similiar effect also exists at the PHENIX central arm 
edges, here due to the lack of balancing fragments. Both issues can be effectively addressed 
by a nonflat weighting that smoothly dampens large angle fragments. Furthermore, the energy 
flow variable as proposed by Sterman et al. 1 3] suggests that angular convolution of the event pj 
with a continuous distribution can provide an effective description of QCD processes. All these 
reasions pointed us towards a seedless jet definition as the output of a filter, e.g. by using the 
Gaussian weighting 

P f\v,<t>)= r r di d<p' p T {i^)e^>-^-* ,f ^ (i) 

*Ji]q J —n 

and reconstructing the resulting local maxima in pf 1 as jets. The Gaussian filter kernel has the 
property of no maximum creation, and the jet definition therefore becomes trivially collinear 
and infrared safe. A description of its efficient implementation and comparison of the behavior 
against other jet definitions for p+p collisions is given in J4j]. 

Historically, there have been similar attempts to define jets using a filter. The British-French- 
Scandinavian collaboration used a jet definition via convolution [5], which notably involves a 
Gaussian kernel with cr = 0.5 and predates the Snowmass accord on the cone algorithm 16]. 
However, little subsequent study and application of the Gaussian filter was made. 

A direct application of jet reconstruction algorithms to heavy ion collisions is known to give 
rise to false apparent jet production rates |7|]. Unlike at LHC, the collision energy at RHIC 
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Figure 1 : PHENIX Run-5 p+p raw jet spectrum using cr = 0.3 Gaussian filter 



results in a significantly lower jet production rate, e.g. for Cu+Cu at y/s^N - 200 GeV and 
p T = 10 GeV, (2n)- l N-^dNI(pTdp T dy) ~ 10~ 6 (GeV/cr 2 (cf.Fig.H Jet production at these 
rates is easily submerged under background fluctuations. The direct fake rejection is an effective 
techique to filter out fake jets, since it provides jet-by-jet information about its likelihood to be 
the result of a real hard scattering. 

This leads us to develop a fake rejection strategy that can achieve a stronger level of fake jet 
suppression than what has been proposed for the LHC (e.g. |8J]), while preserving a fast rise of 
efficiency to unity. We define a fake rejection discriminant 

g<rM)= X p y-«*-vf^-*f)/^, ( 2 ) 

/efragment 

where (77, <f>) is approximately the jet axis, but we allow it to shift, since jets may not be perfectly 
centered due to background fluctuation. The discriminant size is chosen to be cr^ =0.1. 

Fig. Q] shows the reconstructed p+p jet spectrum. Subsequent unfolding using the geant 
derived detector response shows that the spectrum is consistent with [9]. The loss of n and K° L 
shifts the reconstructed energy scale pj C below the true energy scale. The fact that the pj reaches 
out to 40 GeV/c demonstrates that for the study of jet physics, the fast event readout of PHENIX 
easily offsets its small aperture. 

Fig.|2]shows the reconstructed Cu+Cu dijet A(p distribution for different fake rejection thresh- 
olds and symmetric dijets with 7.5 GeV/c < p r j C < 1 1.5 GeV/c. The effect of the fake rejection 
to suppress fake jets that has a random orientation against jets produced by hard scattering is 
clearly visible. For > 4.9 (GeV/c) 2 , the pedestal translates into a fake rate about 10 times 
the jet yield. The effect of the fake rejection saturates at approximately > 17.8 (GeV/c) 2 , 
where the fake contribution is constrained to below 10%. 

Fig. [3] shows the reconstructed Cu+Cu jet spectrum for different centralities and ga-^ > 
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Figure 2: Run-5 Cu+Cu A<f> distribution for symmetric dijets with 7.5GeV/c < Pj C < 11.5GeV/c and different fake 
rejection levels 
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Figure 3: PHENIX Run-5 fake rejected Cu+Cu raw jet spectrum using cr = 0.3 Gaussian filter, > 17.8 (GeV/c) 2 
fake rejection 
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Figure 4: Run-5 Cu+Cu A(p distribution for symmetric dijets with 7.5GeV/c < tfp < 11.5GeV/c, > 
17.8 (GeV/c) 2 fake rejection, and different centralities 

17.8 (GeV/c) 2 . After the application of fake rejection, a consistent power-law shape across all 
centralities is evident. 

Fig.|4]shows the reconstructed Cu+Cu dijet A(p distribution with > 17.8 (GeV/c) 2 and 
for symmetric dijets with 7.5 GeV/c < < 1 1.5 GeV/c in different centralities. Also here, the 
application of fake rejection results in a A(f> distribution that is consistent across all centralities 
(within statistical errors). 

In summary, we showed that jet physics can be effectively studied using PHENIX and the 
Gaussian filter-based jet reconstruction algorithm we proposed. The proper rejection of fake jets 
is an important aspect of applying jet reconstruction to heavy ion collisions. We have shown our 
first results in p+p and Cu+Cu collisions. Further studies, including the measurement of Raa and 
the fragmentation function, are underway. 
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